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We discuss the Auxiliary Field Quantum Monte Carlo (AFQMC) method applied to 
dilute neutron matter at finite temperatures. We formulate the discrete Hubbard- 
Stratonovich transformation for the interaction with finite effective range which is free 
^ ' from the sign problem. The AFQMC results are compared with those obtained from 

CO , exact diagonalization for a toy model. Preliminary calculations of energy and chemical 

OO 1 potential as a function of temperature are presented. 
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1. Introduction 

Properties of dilute and degenerate Fermi gases with large scattering lengths are of 
\ particular interest in the context of cold fermionic atoms and dilute neutron matter. 
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In the case of trapped fermionic atoms, due to very low density, the scattering 
length is the only parameter which determines an atom-atom interaction. Taking 
advantage of Feshbach resonances it can be tuned, using external magnetic field, 
■ to the values which greatly exceed the average interparticle distance. This limit 

which is dubbed the unitary regime exhibits universal properties and has been the 
subject of intensive theoretical effort in the last couple of years (see^and references 
therein). Similarly, neutron matter at densities corresponding to kp < 0.6 fm"-'^ is 
a dilute system of fermions in a sense that the average distance between particles 
is much larger than the range of the neutron-neutron interaction. It is superfluid 
and its physics is determined by two parameters: scattering length as and effective 
range t^q in ^Sq neutron-neutron channel. The influence of other channels as well 
as of three-body forces is marginal and can be neglected in this density rangj2l. 
The properties of dilute neutron matter are crucial for understanding the structure 
and thermal evolution of neutron star crusti^. They can also provide constraints for 
parameters of nuclear energy density functional. 

Despite of these similarities the important qualitative difference between dilute 
neutron matter and cold atomic gases arises from the fact that in the former case 
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the nonzero effective range cannot be neglected. Tlierefore contrary to cold atomic 
gases, physics of dilute neutron matter cannot be captured by simple contact, delta- 
like interaction (except for very low densities kp < 0.06fm~^). 

Since Fermi gas with long scattering length is an example of strongly correlated 
system, non-perturbative methods based on Monte Carlo approaches are often used 
to determine its propert ies . There exist various calculations at zero temperature 
| 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | jjgj.g present selected aspects of AFQMC approach 

at finite temperatures, which is able to generate fully non-perturbative results for 
dilute neutron matter without sign problem. 



2. Theoretical approach 

Neutron-neutron interaction in dilute neutron matter can be described by small 
number of parameters characterizing the low energy physics of two-body collisions. 
In this limit the main contribution to the scattering amplitude comes from s-wave 
scattering and is determined by two parameters: scattering length 05 and effective 
range r^Q i.e.: 

/(P) = w ^1^ (1) 

-1/05 + ~ W 

where as ~ — 18.5/m and Tgg- « 2.8 fm. 

To capture low energy physics of dilute neutron matter we consider the Hamil- 
tonian of the form: 

A— i 

+1 f d'rd'r'v{r-r')^p{{r)i^i{V')i:y{r')4^,{V), (2) 

where the field operators obey the fermionic anticommutation relations 
{ip\{'r*),ip\'{l^')} = dxyS{l^ — ~r") and A denotes the spin degree of freedom. 
Throughout this work we shall use units in which Planck's and Boltzmann's con- 
stants are equal to one. The interaction has the form: 

{6g, ~r' — ~7^' ~ 
g, y - T^' e , (3) 
0, otherwise. 

where A4 = {(&, 0, 0), (-6, 0, 0), (0, b, 0), (0, ~h, 0), (0, 0, 6), (0, 0, -h)}. 
This choice has been motivated, as will become clear later, by the special form 
of discrete Hubbard-Stratonovich (H-S) transformation. It imposes limitations on 
the density range of neutron matter which can be described by this particular 
interaction. Clearly the interaction depends on two parameters: g and h. Solving 
the equation for T-matrix elements: 

T^^, = V-*^, + f V^->G^ ^T->^ (A) 

^ p p' " p p' ^ J ^2Tr'^3 p k p' k k~p" ^ ' 
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(5) 



where G is the free particle propagator and using the relation: 

- — T^L, « - — + ir^fjp^ -ip + Oip") 
m P P as 2 '^'^^ 

we can adjust parameters g and b to reproduce scattering length 05 and effective 
range r^g of the neutron-neutron interaction. Still however the regularization pro- 
cedure is required, which can be determined by introducing the momentum cut-off 
Pcut- This prescription sets to zero all two-body matrix elements, if the relative 
momentum of two particles exceeds a given momentum cut-off. 

In order to compute thermal properties of dilute neutron matter we place the 
system on a 3D spatial lattice. The lattice spacing b and size L — Ngb introduce 
natural ultraviolet (UV) and infrared (IR) momentum cut-offs given by Pcut — T^/b 
and po = 27r/L, respectively. The momentum space has the shape of a cubic lattice, 
with size 27r/& and spacing 27r/L. To simplify the analysis, however, we place a 
spherically symmetric UV cut-off, including only momenta satisfying p < Pcut- 

To evaluate expectation values of observables we have used a path integral 
representation of partition function: 

=Tr|nexphT(i7-M7V)]|, (6) 



0(/3,/i) 



1 



Tr <^ 6 exp[-T(H 

{ k=l 



(7) 



where (3 = 1/T — N^t, T is the temperature and /i - the chemical potential, N 
- the particle number operator and O is a quantity of in terest. The propagator 
e"'^^^^^^-' is decomposed using the second order expansiorP^H^ 



exp[— T(i? — ^iV)] « exp 



r(f - ^lN) 



exp(— ry) exp 



r(f - ^lN) 



(8) 



where T is the kinetic energy operator. To express the many body operator e~^^ 
by a sum of one bod y op erators we have extended the discrete H-S transformation 
introduced by Hirscln^. Namely, one may rewrite the interaction in the form: 



{fi^{~r)n^{~r') + n^{~r)fii{~r') + ni{~r)n^{~r') + 



(9) 



where n>,(~r*) = Note that in the coordinate representation each 

element of the above sum can take the following values: 



(10) 



'0, iffni(y) 




-f rit(V'') 




) = 


0, iff n|(~F*) 








) = 1 


< g/2, iffn|(y) 








) = 2 


3g/2,iff nT(y) 


+ ni{~r) 




+ 


) = 3 


, 6g/2, iff n|( r ) 








) = 4 
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where nx{/r) denotes eigenvalue of density operator. This allows us to introduce 
the H-S transformation of the form: 

g— rV' ^ n - ^ e'^' J- )+"t( )] (11) 

where values of Ui and k have to fulfil the equations: 



1 y^'-- p<y^ ^ 1 

1 y^fc „2(T. 

1 Y^fc 
fe Z^ i=l 
1 v^fc 
fe Z^ i=l 



e ' = e 
fe / ..^-^ = e 



(12) 



This set of equations has a solution if fc ^ 7 and 5 < (attraction) . 
The partition function can be expressed in the form: 

^(A/^)^/!! n 2?^(^,^',T,)TrW(M), 

Z^({a}) = T, cxp{-r[/i({a}) - ^]}, (13) 

where Tt- defines the time-ordered product and h{{ij}) denotes the one-body Hamil- 
tonian. Consequently the expectation value of an arbitrary operator takes the 

J Z{(i,fi) TiUiia}) 

The many-fermion problem has thus been reduced to the typical AFQMC problem. 
Consequently the standard Metropolis algorithm can be applied, using Tr U{{a}) 
as a measurelSI. 

To prove that the measure is positive it is sufficient to note that the operator 
U{{a}) is invariant with respect to time reversal. This is due to the special form 
of expression ((TT|) and the fact that the sigma field is real. Therefore the single 
particle spectrum of U{{a}) consists of pairs which are conjugate to each other. 
Clearly Tr ^Y({cr}) = det(l +Z^({o-})) JJ^ \{l + u,)\'^ > 0. In particular if the single 
particle basis is chosen in such a way that the time reversed partners correspond 
to spin-up and spin-down particles then the matrix representation of U takes the 
form: 

where W|| — tin* . Matrix elements between states with different spin projections 
are equal to zero. Finally one gets 

TiU^ det[l +W] = det[l +^|]det[H-Wf^] = |det[l ^ 0. (16) 

Hence the lack of the sign problem is a direct consequence of the time-reversal 
symmetry which is preserved by each term in the sum (jlip . 
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Table 1. Comparison of AFQMC results with 
exact diagonalization in a restricted space 
(see text for details). T/ep = 0.58, n/ep = 
0.41. 

Observable AFQMC result Exact result 
(1000 samples) 



3. Numerical results 

In order to verify the correctness of our code we have performed several tests. As a 
first check the thermodynamics of a free gas was reproduced for g — Q. 

To check our results for an interacting system we have diagonalized the Hamil- 
tonian (2) exactly, restricting the single-particle Hilbert space to the lowest 7 states. 
It corresponds to the Fock space of 2^'* states. The test provided an estimate for 
the size of the imaginary time step. Comparison of AFQMC results with the exact 
ones for selected observables at fixed temperature is given in the table [T] 

To perform calculations for a neutron matter we have used the box with iV^ = 8 
and lattice constant h ~ 3.21fm. The neutron mass was fixed at 939 MeV. The 
chemical potential was chosen in such a way to keep the total number of particles 
equal to 55. This corresponds to the density of neutron matter p = 0.02po where 
Po is the saturation density of nuclear matter. For this density we have varied 
temperatures from O-OQep (0.26 MeV) to l.Oep (3.8 MeV) where ep is the Fermi 
energy. For the lowest temperature we have used 2360 imaginary time steps while 
for the highest temperature only 216. In all runs the single-particle occupation 
probabilities for the highest energy states were below 0.01 at all temperatures. At 
low temperatures the Singular Value Decomposition technique was used to avoid 
instabilities of the algorithm. 

Figured] presents an example of the results for the energy and chemical potential 
as a function of temperature. The energy and chemical potential versus temperature 
for the free Fermi gas at the same particle density has also been plotted (dotted 
line). The free Fermi gas results were shifted by constant values —0.54 and —0.60 
for the energy and the chemical potential, respectively. 

For our lowest temperature 0.06ep (0.26 MeV) we obtained the energy (relative 
to the energy of the free Fermi gas Effg) E/Effg = 0.48 ± 0.02 {E/N = 1.28 
MeV). This value is in good agreement with recent fixed- node quantum Monte 
Carlo calculationJl^. The chemical potential at this temperature was calculated to 
be n/eF = 0.40 (1.77 MeV). 



(T> 
{V) 

(N) 



0.32 
-0.088 
0.23 
5.58 



0.31 
-0.086 
0.23 
5.49 
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Fig. 1. Energy and chemical potential as a function of temperature for neutron matter at density 
p = 0.02po. Dotted line denotes energy and chemical potential for free Fermi (Free Fermi gas 
results were shifted by constant values- see text for details.). 



4. Conclusions 

In this work we have introduced the simple effective interaction which is capa- 
ble to describe properties of dilute neutron matter. For this interaction we have 
constructed Hubbard-Stratonovich transformation which yield to positive measure 
and allows to perform Monte Carlo calculations without fermionic sign problem. We 
have performed the fully non-perturbative calculations for dilute neutron matter of 
density p = 0.02po at finite temperatures. 
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